NPS71-87-001 (REVISED) 

NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




NOTES ON SEARCH, DETECTION 
AND LOCALIZATION MODELING 



R« N. FORREST 

'/ 



APRIL 1987 

(Revised October 1987) 



Approved for public release; distribution unlimited. 

Prepared for: 

Naval Postgraduate School 
Monterey, CA 93943-5000 



FedDocs 
D 208.14/2 
NPS-7 1-87-00 1/987 




fe-c^L 

Vf)2 

-? - ^ 0 ' ■ ? 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Rear Admiral R. C. Austin 
Superintendent 



K. T. Marshall 
Acting Provost 



Reproduction of all or part of this report is authorized. 



UNCLAbblHLD 



DUDLEY KNOX LIBRARY 



SECURITY CLASSIFICATION OF THIS PAGE 



REPORT DOCUMENTATION PAGE 


la REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED 


lb RESTRICTIVE MARKINGS 


2a. SECURITY CLASSIFICATION AUTHORITY 


3 DISTRIBUTION /AVAILABILITY OF REPORT 

Approved for public release; distribution 
uni imi ted. 


2b DECLASSIFICATION /DOWNGRADING SCHEDULE 


4. PERFORMING ORGANIZATION REPORT NUMBER(S) 

NPS71-87-001 (REVISED) 


5 MONITORING ORGANIZATION REPORT NUMBER(S) 



6a NAME OF PERFORMING ORGANIZATION 

Naval Postgraduate School 



6b OFFICE SYMBOL 
(If applicable) 

Code 71 



7a NAME OF MONITORING ORGANIZATION 



6c. ADDRESS {City, State, and ZIP Code) 

Monterey, CA 93943-5000 



7b ADDRESS (Ofy. State, and ZIP Code) 



8a NAME OF FUNDING / SPONSORING 
ORGANIZATION 

Naval Postgraduate School 



Bb OFFICE SYMBOL 
(If applicable) 



9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 



8c. ADDRESS (Oty, State, and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 



Monterey, CA 93943-5000 


PROGRAM 
ELEMENT NO 


PROJECT 

NO 


TASK 

NO 


WORK UNIT 
ACCESSION NO 


11 TITLE (Include Security Classification) 

NOTES ON SEARCH, DETECTION AND LOCALIZATION 


MODELING 









12 PERSONAL AUTHOR(S) 



Forrest, R. N< 



13a TYPE OF REPORT 


13b TIME COVERED 


14 DATE OF REPORT {Year, Month, Day) 


IS PAGE COUNT 


Technical 


FROM TO 


1987 Oct 


75 



16 SUPPLEMENTARY NOTATION 



17 


COSATI CODES 


18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 


FIELD 


GROUP 


SUB-GROUP 


search modeling, detection modeling, localization modeling 

















19 ABSTRACT (Continue on reverse if necessary and identify by block number) 

The report is an extended outline of some of the course material used in 
courses on search, detection and localization model ing„ 



20 DISTRIBUTION /AVAILABILITY OF ABSTRACT 

OO UNCLASSIFIED/UNLIMITED □ SAME AS RPT □ DTIC USERS 



21 ABSTRACT SECURITY CLASSIFICATION 

Unclassified 



22a NAME OF RESPONSIBLE INDIVIDUAL 

R. N. Forrest 



22b TELEPHONE (Inc/ode Area Code) 

(408)646-2653 



22c OFFICE SYMBOL 

Code 55Fo 



DDFORM 1473, 84 MAR 



83 APR edition may be used until exhausted 
All other editions are obsolete 



security classification of THIS PAGE 

O U.S. Government Printing OHice 1986—606-243 



Preface 



This report is a collection of material that has been used 
in courses on search, detection and localization modeling. Its 
organization follows to some extent material by S. M. Pollock in 
Selected Methods and Models in Military Operations Research which 
is listed in the report bibliography. The report is not intended 
to be a text on these subjects. In particular, in some areas it 
does not provide the depth of coverage that is found in the book 
Search and Detection by Alan R. Washburn which is cited in this 
report as Reference 21. 

Typographical and other errors that have been found in the 
original version of this report have been corrected and several 
sections have been modified or expanded in this revision. 
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I. Detection Models and Signal Detection Theory 

Signal detection theory is the basis for analyzing the 
detection models that are described in this report. In signal 
detection theory, the decision making portion of a detection 
system is called the receiver and a detection experiment is the 
observation by a receiver of input data generated during some 
time interval. The data that is related to a target is called 
signal. The data that is not related to the target is called 
noise. In general, the target data is associated with a 
localization region that in some cases is called a resolution 
cell. When a detection experiment is performed, either the event 
Hq = {the receiver input is noise) or its complement 

= (the receiver input is signal and noise) will occur. In 
the first detection models that are described here, after 
analysis of the input data by a receiver, either the event 
Dq =(the receiver decides the input is noise) or its complement 
D]^ = (the receiver decides the input is signal and noise) will 
occur. Detection models for which Dl is the complement of Dq 
are called binary detection models or forced choice detection 
models. Four events which are important in binary detection 
model are indicated in the Venn diagram of Figure 1. 

The Venn diagram emphasizes a decision problem that is 
associated with a receiver that can be modeled using a binary 
detection model. The problem is this: Under what conditions 

should the event occur? That is, under what conditions 

should a receiver decide that the input data indicates a target 
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was present in an observed localization region either at the time 
or prior to the time of the observation? 



Ho Hi 



Dq Hq 


Dq n Hi 


Di n Hq 


Di n Hi 



Figure 1. Four events of importance in binary detection models. 

In the detection model descriptions that follow, the 
following notation and terminology are used: pf = P(Di|Ho), the 

probability of Di given Hq/ is called the false alarm 
probability; p^j = P(Di|Hi), the probability of Di given Hi, 
is called the detection probability and P = P(Hi) , the 
probability of Hi, is called the prior probability. It is the 
probability that a target will be in the localization region when 
the input data is generated. 

In the detection models, the input to a receiver consists of 
a sequence yi»***ym that is a realization of a discrete 
parameter stochastic process whose random variables are elements 
of a continuous parameter stochastic process. Although the 
receiver input in some cases can be considered to be described by 
a continuous parameter stochastic process, because of the finite 
amount of information (unique data) contained in a bounded 



2 



sequence of finite length, a discrete parameter stochastic 
process is sufficient to represent the receiver input in these 
cases. This is established formally by the stochastic sampling 
theorem. 

In the models, for an observation, the input stochastic 
process has the following characteristics: It is a noise process 

when there is no target data and it is a noise process plus a 
signal process when there is target data. In the first detection 
model that is described in Section III, the signal process is a 
deterministic process. In the second and third detection models 
that are described there, the signal process is a random process. 
To specify the noise process or a random signal process, one 
needs only to specify the joint distribution of the finite 
sequence of random variables that determine the process. If the 
signal process is a deterministic process, the signal values are 
a sequence of values that are assumed to be known before an 
observation is performed. Consequently, in this case, to specify 
the process, one needs only to specify these values. 
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II. Decision Criteria 

To simplify the discussion of decision criteria and 
decision rules, a receiver's input will be assumed to be 
determined by a single decision random variable Y. In this 
case, the input process in determined by the conditional 
distribution function Fy(y|Ho) when the input is noise alone 
and by the conditional distribution function FyCylHj^) when the 
input is signal plus noise. 

The condition that a receiver's input is required to 
satisfy in order that the event will occur can be specified 

in temns of a decision rule. For the assumed case, a decision 
rule is a rule which determines for every observable value of Y 
the decision that the receiver is to make. The decision rule can 
be considered to be a function <P(y) which relates each 
observable value y to one or the other of the following two 
commands: do = "Decide that the receiver input was noise." and 

dj^ = "Decide the input was signal and noise.". Choosing a 
decision rule <fi(y) defines a set D of observable values of Y 
such that the event D^ = { Y e D ) . 

The problem which was considered in Section I can now be 
restated in the following way: What criterion should be adopted 

in order to determine a decision rule or, equivalently, its 
corresponding set n ? A desirable characteristic for a criterion 
is suggested by the following argument: Consider the odds in 

favor of given y is observed. That is, consider 
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P(H 3 ^|y = y)/P(Ho|Y = y) . One might expect that y would be a 
member of the set fl if and only if y made this ratio equal to 
or greater than some value k. But this is equivalent to defining 
n as follows; n = { y : L(y) > K ) where L(y) is the 
likelihood ratio associated with an observed value y and K is 
a constant related to the constant k. This suggests that 
choosing an optimum criterion is equivalent to choosing an 
optimum value for K. Four specific decision criteria are defined 
next in terms of K. For each criterion, fl has the above 
form. But for each criterion the choice of K is different. 

The decision criteria are; 

1. The Neyman-Pearson Criterion; Choose n so that p<j 

is a maximum subject to the constraint that pf < a where a is 
a specified value. For a continuous decision random variable, 
the constant K is chosen so that pf = a. 

2. The Bayes Criterion: Choose n so that the expected 

cost of a receiver's decision is a minimum. For a continuous 
decision random variable, K = [ (cio'^^oo)/ ('^01“'^ll) ] * 

°10 ^ ^00 *^01 ^ *^11 where Cj^j is the cost of Dj^ n Hj . 

3. The Ideal Observer Criterion; Choose n so that the 
probability that the receiver makes an incorrect decision is a 
minimum. K = (1-P)/P for a continuous decision random variable. 

4. The Minimax Criterion; Choose fi when P is unknown so 

that the maximum expected cost of a receiver's decision is a 
minimum. If c^q > *^00 ^01 ^ '^ll' then 
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K = [ (cj^o”<^00)/ ^ continuous decision 
random variable. Here, P* is the value of the prior probability 
P that would make the expected cost of a receiver's decision a 
maximum if P were known and the Bayes Criterion were used. 

If a model which specifies the conditional distributions 
Fy(y|Ho) and FyCylHj^) and a decision rule are adopted, then 
the value of pf and the value of p^j are determined. This 
pair of values (Pf»Pd) called a receiver operating point. 

If the decision rule results from using a likelihood ratio 
criterion such as one of the four listed above, then it will 
involve the parameter K since fl = {y: L(y) > K). And, for a 
given value of K, since n uniquely determines the pair 
(Pf,Pd) , a single operating point results. By varying K, a set 
of operating points can be generated which determines a receiver 
operating characteristic curve or ROC curve. Different ROC 
curves can be produced by changing either one or both of the 
conditional distributions which implies changing either the 
signal process or the noise process. 

A decision rule which results from using a likelihood ratio 
criterion in a model in which the input process is determined by 
a set of m random variables can be expressed in terms of a set 
n as follows: a = { (yi,...,yn) : L(Yi» • • • ^Ym) ^ K ) where K 
is specified in the same way as in the case in which m = 1. 
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III. Three Binary Detection Models 

Three detection models are examined in this section. 

For the first two detection models, the input stochastic process 
for an observation is defined by a time sequence of continuous 
random variables. The random variables represent a sample from a 
continuous parameter stochastic process which is sampled at times 
such that the random variables are independent. For the third 
detection model, the input stochastic process is a counting 
process and it is defined by a single discrete random variable 
that is equal to the number of events that are counted during the 
observation. 

Model I: In the first detection model, a sampled noise 

value is a value of a normally distributed random variable with 
mean zero and with known variance a’ and a sampled signal value 
is a known value of a deterministic variable. Thus, the input 
process corresponding to an observation consists of some number 
m of independent normal random variables Y 2 ,...,Vj„ each with 
variance a’ . And, for i = 1,2,3, ••• , m, when a signal is not 
present the mean of is zero and when a signal is present the 

mean is s^. The result of using a likelihood ratio decision 
rule in the model can be expressed in terms of a random variable 
Z. This random variable is called a crosscorrelation statistic 
and it is defined by Z = S sj^*Yj^ where the sum index 1 = 1,2, 
,m here and in the remainder of this section. However, it 
is more convenient to express the result in terms of a statistic 
V which is defined by V = In terms of this random 
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variable, the conditional probabilities pf and p^j are given 
by: Pf = 1 - #(v*) and Pd = 1 ~ *(v* - d^) where ♦ 
symbolizes the standard normal ctimulative distribution function, 
the threshold v* = (l/a 2 )*(a^ln K + (1/2) E s^^) and 
d = Z s,\/a^ . The parameter d is called the detection index. 

Often, the input stochastic process represents a quantity 
whose square is proportional to power. In such a case, the 
average receiver input power is the random variable Z Y^/m. If 
a signal is not present, the expected average receiver input 
power is N = Z aV® = where N is called the noise. The 
average receiver input signal power is S = Z s^/m where S is 
called the signal. In these terms, d = m* (S/N) where S/N is 
called the signal-to-noise ratio. 

If a receiver's input can be considered to be a time 
sequence of continuous voltage values such as in the case of a 
sonar receiver, in some cases a frequency representation can be 
used that involves the concept of receiver bandwidth. In these 
cases, the noise process is assumed to be such that m = t/£t 
where t is the integration time (the duration of an 
observation) and St is the time between samples with 
<St = 1/[2(BW)] where BW is the bandwidth and St is 
determined by the sampling theorem. This implies that the 
detection index can be written as d = 2t« (BW) • (S/N) . By 
defining Nq as the power spectral density where Nq = N/BW^ the 
detection index can also be written as d = 2t* (S/Nq) . In 
Reference 2, the conditions required by the first model are 
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called Case I and in the following sections the first model is 
called the Case I model. A receiver that processes data such 
that it would implement a likelihood ratio decision rule under 
the conditions of the first model is called a crosscorrelation 
detector. If the description of the input noise is adequate, a 
Case I model can be used to obtain an estimate of an upper bound 
on a detection system's performance, since all the information 
necessary to define the signal is assumed to be known. 

Model II: In the second detection model, a sampled noise 

value is a value of an independent normal random variable with 
mean zero and with variance a' . And a sampled signal value is 
an independent random variable with mean zero but with variance 

Thus, the input process corresponding to an observation 
consists of some number m of independent normal random 
variables Y3^,...,Yjj, each with mean zero and each with variance 
when a signal is not present and with variance a' + when 
a signal is present. The result of applying a likelihood ratio 
decision rule in this model can be expressed in terms of a 
statistic X which is defined by X = E Y^. When a signal is 
not present, the statistic X/N has a chi-square distribution 
with m degrees of freedom. When a signal is present, the 
statistic X/(N+S) has a chi-square distribution with m 
degrees of freedom. So, in terms of these statistics, the 
conditional probabilities pf and p^j are defined as follows: 
Pf= P(X^ > x*/N) and Pd = P{XJ, > (x*/N) -[l/Cl+S/N)]) 
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where is a chi-square random variable with m degrees of 

freedom, x* is a number which is determined by the decision mile 
and S/N is the signal-to-noise ratio. A receiver that would 
implement a likelihood ratio decision rule under the conditions 
of the second model is called an energy detector (or square law 
detector) . 

The mean of a chi-square random variable with m degrees of 
freedom is m and the variance is 2m. By the central limit 
theorem, as the number of degrees of freedom of a chi-square 
random variable becomes large it can be approximated by a normal 
random variable with the same mean and variance. Hence, for a 
sufficiently large sample size m or equivalently, for a 
sufficiently large time bandwidth product t* (BW) , the 
conditional probabilities pf and p^j can be approximated by 
Pf = 1 - *(v*) and Pd = 1 - ♦{ [1/(1+S/N) ] • (V* - d^) ) where the 
threshold value v* = (x* - ma* )/( (2m) ] and d = t*(BW)(S/N)^ 
if the concept of bandwidth is applicable. If, in addition, the 
noise N is significantly larger than the signal S, that is, if 
t»(BW) » 1 and S/N « 1, then pf and p^j can be 
approximated by: pf = 1 - *(v*) and Pd = 1 ~ #(v* - d^) . 

In Reference 2, the conditions required for this approximation 
are called Case II and in the following sections the limiting 
form of the second model is called the Case II model. In 
Reference 2, an example is given that uses the Case II model to 
estimate the performance of a passive sonar system. 
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Model III; In the third detection model, a sampled noise 
value and a sampled signal value are values of independent random 
variables that are determined by independent Poisson processes 
observed for a time interval t. The noise process is 
characterized by a counting rate a, the signal process is 
characterized by a counting rate Og and the noise and signal 
processes are additive. Consequently, when the input is noise 
alone, the input is the value of a random variable Y that has a 
Poisson distribution with parameter a*t. And, when a target is 
not present, the input is a value of a random variable that has a 
Poisson distribution with parameter (a + ag)*t. For a 
likelihood ratio decision rule, pj = 1 - P(y*;a»t) and 
P(j = 1 - P[y*;(a + Og) *t] where P(y;0) symbolizes the 
cumulative Poisson distribution function with parameter 0 and 
y* is a threshold value of y determined by the decision rule. 

In terms of the above notation, both the mean and variance of a 
Poisson random variable are equal to 0. When 0 is large, the 
cumulative Poisson distribution function can be approximated by 
the cumulative distribution function of a normal random variable 
having the same mean and variance. Using this approximation for 
cases where a*t is large gives pf = 1 - *(v*) and 
Pd = 1 - ♦{[1/(1 + ag/a)^](v* - d^)} where v* = (y* - at)/(at)'^ 
and d = a*t*(ag/a)’. In addition, if the signal counting rate 
is significantly smaller than the noise counting rate, that is, 
if a»t » 1 and Qg/a « 1, then pf and p^j can be 
approximated by: pf = 1 - #(v*) and Pjj = 1 - #(v* - d^) . 



11 



This detection model might be used to describe a receiver whose 
input for an observation is the number of photons counted by a 
radiation detector in situations where a«t, the expected number 
of counts when no signal is present, is of the order of thirty or 
more. In this case, the target might be a wake segment for 
example. 

When a likelihood ratio decision rule is used in the three 
models discussed above, for the first model and under limiting 
conditions for the second and third models, the following result 
is obtained; pf = 1 - #(v*) and = 1 - #(v* - d^) where the 
definition of v* depends on the noise power N for the first 
and second models. For a sonar receiver described by the first 
model, that is, by the Case I model: d = 2t» (BW) (S/N) . For a 

sonar receiver described under the limiting conditions for the 
second model, that is, by the Case II model, d = t • (BW) (S/N) * . 
So, in either a Case I model or a Case II model of a sonar 
receiver, the detection index d is a function of the time 
bandwidth product t»(BW) and signal-to-noise ratio S/N. Since 
sonar equations relate S/N to system, target and environmental 
parameters, a sonar equation can be used to relate S/N to these 
parameters in a model of a sonar receiver. 
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IV. General Detection Models 

The detection models that have been considered to this 
point are based on binary detection theory. After each 
observation, a receiver decides either that the input 
corresponding to the observation was noise or else it decides it 
was signal plus noise. However, in some detection systems this 
decision is delayed. In a computational sense, a model of such a 
detection system is generally complex relative to a binary 
detection model. To illustrate this, consider an active sonar 
system whose receiver includes an operator. Suppose the 
probability that the operator will detect a target echo has been 
determined in a laboratory experiment in which the operator was 
required to decide after each input corresponding to a resolution 
cell that either the input was a target echo (signal) and noise 
or the input was noise alone. In addition, suppose that under 
operational conditions the operator normally delays this 
decision. Then, in general, the probability that the operator 
will decide that the input corresponding to a resolution cell 
that contains a target is a target echo and noise will not be 
equal to the probability of the event in the forced choice 
experiment. And, in addition, the probability that the operator 
will decide the input corresponding to a resolution cell that 
does not contain a target is a target echo and noise will not be 
equal to the probability of this event in the forced choice 
experiment. Consequently, in general, the value of both p<j and 
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Pf for an operational environment will be different than that 
for the laboratory environment. 

One model that has been proposed to deal with this kind of 
situation defines the event that a receiver decides that the 
input corresponding to a resolution cell is signal and noise to 
be equivalent to the event that out of n consecutive 
observations at least k of them would result in the decision 
that the input was signal and noise in a forced choice 
experiment. The model is said to be based on an k-out-of-n 
detection criterion. With this criterion, the probability that a 
target will be first detected on the observation can be 
found as follows: Determine the 2^ sequences of forced choice 

responses that could result for a sequence of j consecutive 
observations. Next, determine the probability of occurrence for 
each sequence that first satisfies the k-out-of-n detection 
criterion on the observation. The probability of first 

detection on the observation is equal to the sum of these 
probabilities. The cumulative probability of detection at the 
jth observation is the sum of the probabilities of first 
detection on the i^^ observation for i=l,2, ••• ,j. 
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V. Signal-to-Noise Ratio Detection Models 

In some radar and sonar detection models, for a 
specified value of pf, a minimum acceptable value of p^j is 
defined. This minimum acceptable value of p^j and the specified 
value of Pf define what can be called a minimum acceptable 
signal-to-noise ratio (S/N)jj if p^ is a nondecreasing 
function of S/N. In some sonar detection models, this 
signal-to-noise ratio in decibels is called the detection 
threshold DT. In symbols, DT = 10 log(S/N)]g. If the minimum 
acceptable value of p<j is .5, then DT is usually called the 
recognition differential RD. The difference between the 
signal-to-noise ratio in decibels and RD (or DT) is called the 
signal excess SE. In symbols, SE = 10 log(S/N) - RD. 

One interpretation of signal excess is that for a 
localization region containing a target detection occurs with 
probability one if SE > 0 and with probability zero if SE < 0. 
This interpretation provides the basis for defining detection in 
the three encounter detection models that are discussed in 
Section VII. A more consistent interpretation of signal excess 
is: If SE > 0, then the probability of detection p^j is 

greater than or equal to the minimum acceptable value (.5 if 
recognition differential RD is used to define signal excess) . 
For cases where p^j increases rapidly with signal excess in the 
neighborhood of zero signal excess, the two interpretations are 
essentially equivalent. For a discussion of this point as well 
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as a discussion of an operational case In which receiver 
decisions are delayed, see Reference 3. 

Signal excess (slgnal-to-nolse ratio) detection models 
provide a basis for general detection models. In particular, 
models that describe nonstationary noise and signal processes and 
randomly changing decision rules. This Is Illustrated by the 
models described In Section VII. In addition, signal excess 
models provide a basis for delayed receiver decision models. 

This Is Illustrated by the active sonar detection models In both 
Reference 4 and Reference 5 that are based on a k-out-of-n 
detection criterion. In all of these models, the slgnal-to-nolse 
ratio and the recognition differential are random variables. 

Using X(t) to represent a random variable corresponding to 
an Index time t and a subscript to Identify the random variable 
In such models, for a passive sonar receiver, the signal -to-nolse 
ratio In decibels associated with a decision at the Index time 
Is; XsL(t) " X<i>L(t) “ [Xnl(^) " Xj 3 i(t)). In this expression, 

SL represents source level, TL represents transmission loss, 

NL represents noise level and DI represents directivity Index. 
Since signal excess SE Is defined to be the difference In 
decibels between the signal-to-noise ratio and the recognition 
differential (or detection threshold) , it too is a random 
variable and, for any decision time t, one can .write: 

(1) XsE(t) = XsL(t) - XTL(t) - [XNL(t) - Xqi (t) ] - XRo(t). 

The distributions of the random variables on the right side of 
Equation 1 determine the distribution of the signal excess. In 



16 



the passive sonar detection model described in Reference 6, 

' XRD(t) and, in effect, normally 

distributed random variables while XipL(t) is a uniformly 
distributed random variable. In the three signal excess models 
that are described in Section VII, all of the random variables in 
Equation 1 are normally distributed. 

It is sometimes convenient to write Equation 1 as follows: 

(2) XsE(t) = SE(t) + X(t) . 

In Equation 2, SE(t) is the expected value of the signal excess 
determined by the following expected value equation; 

(3) SE(t) = SL(t) - TL(t) - [NL(t) - DI(t)] - RD(t) 

where each term on the right represents the expected value of the 
indicated random variable and X(t) is a random variable that 
determines the stochastic character of the signal excess. Since 
SE(t) is the mean of XsE(t)» Equation 2, the mean of X(t) 
is equal to zero and the standard deviation of X(t) is equal to 
the standard deviation of XgE(t) . If o represents the 
standard deviation of Xse(^) random variables on the 

right side of Equation 1 are statistically independent, then 

“ ‘^SL ^TL ‘^DI ‘^RD* This relation has been used to 

determine a standard deviation for the signal excess in 
operational models. 
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VI. General Encounter Models 



A basic problem associated with search modeling is that 
of determining the probability that a target will be detected by 
a detection system during an encounter with one or more detection 
systems. In the encounter models that are considered in this 
report, during a search, observations are made of a series of 
localization regions. For a localization region that contains a 
target, P(j is the probability of the event {D^ n H^}, that is 
the event detection. For a localization region that does not 
contain a target, pf is the probability of the event (D^ n Hq), 
that is, the event false alarm. In addition, in these models, 
the time to resolve a false alarm is ignored. However, Pjj and 
Pf are assumed to be determined by some criterion such that Pf 
is less than one. 

Using the order number of a decision rather than its time as 
an index relative to detection decisions for localization regions 
that contains a target and a random variable N to represent the 
decision order number at which detection first occurs, the 
probability of detection during an encounter can be written as: 

P(N < n) = P(N < ■) + P(N = m+1) + ••• + P(N = n) or equally as 

P(N < n) = 1 - [1 - P(N < m)]-(l - gjH-i) (1 - gn) where 

g^ = P(N = i|N < i-1) is the probability of the event detection 

at the i^^ decision conditioned on the event no, detection at an 
earlier decision and 1 < m < n. 

The second expression above for P(N < n) is generally of 
greater interest than the first expression, since g^ can 
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usually be more directly related to operational parameters such 
as target range and target environment that determine a target's 
detectability than can P(N = i) . If a target is not present in 
the localization region that corresponds to a decision, then the 
event {D^ n can not occur and the probability of this event 
will be zero. 

With a time rather than the order number to index a decision 
and a random variable T to represent the time index at which 
detection first occurs, P(N < n) becomes P(T < t^) with 
P(T < t„) = 1 - [1 - P(T < tm)]-[l - g(tn,+i)] ••• [1 - g(tn)] 
where now gj^ = P(T = t^lT . 

If g(t£) « 1 for 1 = 1,2, ••• ,n, then, to a first 
approximation, ln[l - g(ti)] = -g(ti) for 1 = 1,2, ••• ,n and 
P(T < tn) = 1 - [1 - P(T < tm)]*exp[-E g(ti)]. This follows 
since P(T < t^) = 1 - [1 - P(T < tm)]*exp[E ln[l - g(ti>] where 
the sum index i = m+1, ••• ,n. A continuous analog to this 
approximation can be used to describe an encounter for which 
g(ti) « 1 for i = m,m+l, ••• ,n and decisions during the 
encounter can be considered to occur continuously. That is, the 
time of an observation corresponding to a decision and the time 
between decisions are both negligible relative to the time of the 
encounter. 

The analog can be developed as follows: First, let fit be 

the time between decisions, then t^ = i*fit and the probability 
of detection P(T < tn) = 1 - [1 - P(T < tni)]*exp[-£ T(ti)*fit] 
where t ( t^) = (1/fit) *g(ti) is a detection rate function 
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(a probability of detection per unit time) and, in terms of fit, 
the probability g(tj^) = P[T = i*fit|T < (i-l)*fit]. 

If T is considered to be a continuous random variable, the 
expression for P(T < t^) above indicates that the sum in the 
exponent should be replaced by an integral whose integrand is a 
continuous function T(t). If T(t) can be determined, then, 
with g(ti) as a guide, the cumulative probability of detection 
P(T < t) can be defined by: 

(4) T(t) = lim { (1/fit) *P(t < T < t+fit|T < t) ) 

where the limit is for fit approaching zero. Equation 4 
implies the differential equation: dp(t)/dt = [1 - p(t)]*r(t) 

where p(t) = P(T < t) . The solution to this differential 
equation is: 



where t is the time index for a decision during an encounter, tm 
is some time during the encounter and tn > tm. A f(t) that is 
based on a visual detection model is described in Reference 7. 

If the detection capability of a detection system is assumed to 
depend on a target's position relative to the detection system 
during an encounter but not to depend on the clock time, then the 
time index of a decision can be a relative index that determines 
the target position that is associated with a decision rather 
than the clock time associated with the decision. 

The above results apply to the case of an encounter between 
a target and a collection of detection systems. However, if the 
detection systems are not collocated, it is generally convenient 



(5) 



P(T < tn) 
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to describe encounters of this kind in terms of encounters 
between the target and the individual detection systems. In 
either case, if the event target detection for a detection system 
is not independent of the event for other detection systems, then 
in order to describe this in an encounter model the correlation 
between the input to the detection system and the inputs to the 
other detection systems must be specified. This has been done in 
some models as follows: First determine the probability of 

detection for each system acting alone. Let Pj^ be the 
probability that the i^^ system detects the target during the 
encounter under this condition. Next, consider two cases: In 

the first case, the random factors that determine detection for a 
system are independent of those that determine detection for the 
remaining systems. In the second case, the random factors that 
determine detection for the systems are completely dependent. In 
the first case, the probability that at least one system detects 
the target is given by: Pj = 1 - (1 - P^) • (1 - P2) *•* (1 ~ Pn) 
where n is the number of detection systems involved. In the 
second case, the probability that none of the systems detect the 
target is given by: 1 - Pq = 1 - Pj, where P^ > Pj^ for 
i = 1,2, ••• ,n since if the m^^ system does not detect the 
target, none of the remaining systems will detect it. The 
probability that at least one system detects the target is given 
by: P = o*P0 + (1 - a) ‘Pj where a determines the degree of 
correlation and 0 < a < 1 . A way to determine a value for a 
is described in Reference 8 . 
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VII. Three Signal Excess Encounter Models 

In the three models described in this section, 
detection is defined in terms of signal excess as it is in 
Section V. Each of the models determines a cumulative 
probability of detection for a target in an encounter with a 
passive sonar system. An observation in the models is indexed by 
time and the index can usually be considered to be the time at 
the end of the observation. During an encounter, observations 
are made of a series of localization regions. For localization 
regions that do not contain a target, the signal observed during 
the observation of the region is zero. The time to resolve a 
false alarm is zero. However, false alarms are not ignored in 
that the value of RD (or DT) is determined by some specified 
false alarm probability. Consequently, although the time to 
resolve a false alarm is zero, the cost associated with a false 
alarm is not zero. In determining signal excess in the models, 
it is convenient to use Equation 2. For each decision in an 
encounter, there is a random variable X(t) defined by Equation 
2 that determines the random character of the signal excess. 

For a sequence of decisions, the set of these random variables 
ordered by their time index constitutes a stochastic process. 

The joint distributions of these random variables determines the 
nature of this process. In the three encounter models described 
in this section, the stochastic process is called a lambda-sigma 
jump process. The jumps in the series occur at times determined 
by a Poisson process with a mean rate lambda. This implies that 
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the time between jumps is a random variable with an exponential 
distribution and that the expected times between jumps t is 
equal to the reciprocal of lambda. The time series that are 
generated by lambda-sigma jump processes are represented by 
Figure 2 below. 




Figure 2. A time series representing a realization of a lambda- 
sigma jump process. On the plot, a in dB equals one unit on 
the vertical axis and r equals one time unit on the horizontal 
axis. 



From Figure 2, note that the observed values of neighboring 
random variables are equal unless a jump has occurred between 
them. When a jump occurs, the first random variable after the 
jump is normally distributed with mean zero and variance a’ and 
it is independent of all the random variables before the jump. 
Conditioned on a jump pattern, this random variable and all the 
random variables between it and the next jump are dependent and 
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the correlation coefficient between any pair is one. That is, if 
the value of the signal excess is known at some time, then all of 
the values between the last jump before that time and the first 
jump after that time are also known. However, since the jumps 
occur randomly, knowing the value of the signal excess with 
certainty at some time does not determine the values of the 
signal excess with certainty at neighboring times. In the 
unconditioned case, the correlation coefficient between the 
random variables X(t) and X(t+r) is equal to 1/e. For this 
reason, r is referred to as a relaxation time. If only the 
transmission loss is considered to be determined by a stochastic 
process, then only the transmission loss will be a random 
quantity and r is the expected time between jumps in the random 
component of the transmission loss. It appears that the use of 
the lambda-sigma jump process is based more on past practice than 
on experimental justification. In this regard, see Reference 9. 
By referring to Equation 1, it can be seen that the lambda-sigma 
jump process is equal to the sum of the stochastic processes that 
determine the random variables on the right side of this 
equation. But, the nature of the processes should depend on the 
encounter that is being modeled. This suggests that an adequate 
description of an encounter using a lambda-sigma jump process 
might not be possible for some encounters. 

In the three encounter models described below, detection is 
defined in terms of signal excess as described in Section IV and 
decisions are indexed by a time that can usually be considered to 
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be the time of the decision. During an encounter, observations 
are made of a series of localization regions. For localization 
regions that do not contain a target, the signal observed during 
an observation of the region is zero. In the models, the time to 
resolve a false alarm is zero. However, false alarms are not 
ignored in that the value of RD (or DT) is determined by some 
specified false alarm probability. Consequently, although the 
time to resolve a false alarm is zero, the cost associated with a 
false alarm is not zero. 

The First Passive Sonar Encounter Detection Model: This 

model describes an encounter in terms of a series of decisions 
with each decision based on the signal excess Xsg(t) at a time 
corresponding to the end of an observation. The observations are 
of equal duration and the integration time that determines the 
recognition differential is equal to the duration of the 
obseirvations. In the model, XsE(t) is determined by a lambda- 
sigma jump process . For an encounter in which SE(t) is 
unimodal and in which the time of the single maximum is prior to 
the end of the encounter, it is shown in Reference 10 that the 
probability p that detection will occur during the encounter 
is given by the following equation: 

(6) p = 1 - [(1 - Pc)/(1 - ^*Pc)]*(l ■ /9*Pi) **• (1 " ^*Pm) 
where ^ = 1 - exp ( -fit/ t) and Pi = *[SE(ti)/a] for i = 1,2, 
••• ,m. Here, fit indicates the duration of an observation and 
♦ indicates the standard normal cumulative distribution function 
as before. The integer c is the index of the decision time 
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for which SE(tc) is greater than or equal to SE(tj^) for any 
time ti and tj < t^ ^ 

As r approaches zero, 0 approaches one and Equation 6 
approaches this form: 

(7) p = 1 - (1 - Pi) ••• (1 - Pxn) . 

In this limit, the signal excess random variables are all 
independent. Note that Equation 7 applies without the condition 
that SE(t) be unimodal. 

As T approaches infinity, 0 approaches zero and 
Equation 6 approaches this form: 

(8) P = Pc- 

In this limit, the correlation coefficient between any pair of 
signal excess random variables is equal to one. Note that 
Equation 8 applies without the condition that SE(t) be 
unimodal. Equation 8 defines a complete dependence encounter 
model . 

The Second Passive Sonar Encounter Detection Model: This 

model is in a sense a third limiting form of the first passive 
sonar encounter detection model. In this limit, the time between 
decisions approaches zero. However, in this limit the 
integration time that determines the recognition differential is 
not equal to fit and it does not approach zero. It is, in 
effect, chosen by the user of the model through the user's choice 
of the value for the recognition differential. For an encounter 
that begins at tl and ends at t2 and for which X 3 £(t) is 
determined by a lambda-sigma jiimp process and SE(t) is 
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unimodal, it is shown in Reference 10 that for this limit. 
Equation 6 has the following form: 



(9) p = 1 - [1 - p(tc)]*exp[-(l/T) . _ p(t) dt] 

where p(t) = ♦[SE(t)/a] and where now t^, is the encounter 
time such that SE(t(,) is greater than equal to SE(t) for any 
other encounter time t and tl < t^ ^ t2 . 

The Third Passive Sonar Encounter Detection Model: This 

model describes an encounter between a target and a passive sonar 
detection system in which detection occurs during an encounter if 
the average value of the square of the continuously observed 
signal-to-noise ratio over a time interval of length u is 
greater than or ecjual to the square of the signal-to-noise ratio 
that determines the recognition differential for an integration 
time equal to u. With R(s) the random signal-to-noise ratio 
at a time s and Rnj(u) the random signal-to-noise ratio that 
determines the random recognition differential for an integration 
time u, detection during an encounter occurs at the first time 
t that the following inequality is satisfied: 



where the time origin is chosen so that t > 0 and where the 
integration time u = t for t < to and u = to for t > to 
where to is a maximum integration time. The random integrand 
in the inequality is related to the random signal excess at the 
time s for an integration time u. The relation is: 




(10) (1/u) _u [R(s)/Ru,(u) ]> ds > 1 




(11) 10 log [R(s)/Rjn(u) ]' = 2[SE(s;u) + X(s) ] 
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where SE(s;u) is the expected value of the signal excess at a 
time s for an integration time u and X(s) is the random 
component of the signal excess at the time s. In the model, 

X(s) is determined by a lambda-sigma jump process and SE(s;u) 
is determined by an expected value sonar equation with a 
recognition differential RD(u) = 10 log rjn(u) . Here, ^^(t) is 
the value of the signal-to-noise ratio that gives a probability 
of detection equal to .5 for an integration time t and a 
specified probability of false alarm pf. With the signal 
detection process described by a Case II signal detection model, 
the detection index necessary to give the required operating 
point (pf,.5) is related to the integration time t and the 
signal-to-noise ratio rjjj(t) by; 

(12) d = u-(BW)trjn(t)]^ 

where BW is the bandwidth of the receiver. For a spectrum 
analyzer, BW would be the bandwidth corresponding to a given 
frequency resolution and d would be the detection index 
required in order to be at the operating point (pf,.5) for a 
signal that was contained within a bandwidth BW. Since d in 
Equation 12 must be the same for t = u and t = to, 

(13) RD(u) = 5 log(to/u) + RD(to) 

where to is the maximum integration time. Then, since 
SE(s;u) - SE(s;to) = RD(to) “ RD(u) , by using Equation 13 and 
Equation 11, Relation 10 becomes: 

ft (1/5) [X(s) + SE(s;to) - 5 log(to)] 

(14) Jt-u 10 ds > 1 
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where as above the time origin is chosen so that t > 0, the 
integration time u = t for t < tQ and u = to for t > to 
and where SE(s;to) is the expected value of the signal excess 
at the time s for a recognition differential determined by an 
integration time to* In an encounter, detection occurs the 
first time that Relation 14 is satisfied. 

As is pointed out in Reference 11, the appeal of the Third 
Passive Sonar Encounter Detection Model relative to the Second 
and First Passive Sonar Encounter Detection Models is that it 
appears to more closely describe the detection process in passive 
sonar detection systems that display their processed data to an 
operator in a continuous manner over a time window of duration 
tQ. However, results reported in Reference 12 indicate that the 
difference between the three models may not be significant in 
some types of encounters. 
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VIII. straight Line Encounters 

Suppose a target's detectability depends on its range 
from a detection system and that the probability of detection is 
effectively zero beyond a range rj^ for any target azimuth. In 
this report, an encounter between the target and the detection 
system is the event that the range between the target and the 
detection system is less than or equal to rj^. In addition, 
suppose rjQ is small enough so that when the target and the 
detection system are having an encounter they can be considered 
to be moving on planes parallel to a tangent plane to the earth's 
surface at some point in their vicinity. If this is the case, 
then while the target and detection system maintain a constant 
course and speed during an encounter, the encounter is called a 
straight line encounter. 

A straight line encounter can be described in terms of a two 
dimensional rectangular coordinate system whose plane is parallel 
to the tangent plane to the earth. If the coordinate system is 
stationary relative to the detection system with the detection 
system located at the origin and is oriented so that the target's 
motion is parallel to the y-axis and is in the positive 
y-direction, then the target's x-coordinate during a straight 
line encounter will be constant. The constant is equal to the 
target's horizontal range at the closest point of approach (CPA) 
on the straight line track on which the target is moving relative 
to the detection system during the encounter. This range is 
called the target's lateral range. 
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A complete straight line encounter is a straight line 
encounter that begins at a range from a detection system that is 
greater than or equal to r^^ and continues past CPA to a range 
from the detection system that is again equal to or greater than 
rj„. Let p(x) be the cumulative probability that a target is 
detected by a detection system in a complete straight line 
encounter in which the target's lateral range is x. Then the 
function p(x) defines what is called a lateral range curve or 
lateral range function. 

Let p be the probability that a target is detected during a 
complete straight line encounter. If the lateral range of a 
target in a straight line encounter is assumed to be a continuous 
random variable X with a uniform distribution with 
for |x| < a/2 and p(x) = 0 for |x| > a/2, then the 
probability that a target will be detected during a complete 
straight line encounter is given by: 



where the limits of integration can be used since the value of 
p(x) is zero for |x| > a/2. Equation 15 suggests a measure of 
a detection system's capability to detect a target in a straight 
line encounter. The measure W is called sweep width and 



(15) 




dx 



(16) 




With this definition, Equation 16 becomes: p = (l/2a) *W. 



31 



IX. Two Intermittent Signal Encounter Models 

In the intermittent signal encounter models that are 
described in this section, an encounter is a complete straight 
line encounter, and during an encounter a target either emits a 
signal (an acoustic transient) or its presence (a visible 
submarine mast) is the cause of a signal at various times. Two 
cases are considered: In the first case, the signals occur 

periodically, the signals are of length 5t and the time between 
the occurrence of signals is r where t > fit. In the second 
case, fit = 0 (the signals are instantaneous) and the signals 
occur at times determined by a Poisson process for which the 
expected time between signals is equal to t. In the model, the 
detectability of a target signal depends on a target's horizontal 
range from a detection system, but on no other factors. If a 
signal occurs while the target is within a range r, it will be 
detected. For a continuous signal, the lateral range function of 
a detection system for a target is: p(x) = 1 for |x| < r and 

p(x) = 0 for |x| > r where the horizontal range r is 
determined by the characteristics of the detection system and the 
target. The geometry for an encounter is shown in Figure 1 
below. 

For intermittent signals, the length of a target's track 
relative to a detection system on which a signai will be detected 
is 2*(r^ - X' ) ^ + wfit where w is the speed of the target 
relative to the detection system. So, a target's exposure time 
during an encounter is (2/w)»(r* - x* ) ^ + fit. 
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For periodic signals, there are two cases. In the first 
case, r > w*(t - 6t)/2. In this case, the signals result in the 
following lateral range function: 
p(x) = 0 for |x| > r 

(17) p(x) = 1 for |x| < (r^ - [w(T - 6t)/2]M*^ 
p(x) = [2/ (w f ) ] • (r^ - x^ ) ^ + fit/r otherwise 

In the second case, r < w*(t - 5t)/2 and the middle equality 
in Equation 17 does not apply. 

For signals that are instantaneous and whose occurrence is 
determined by a Poisson process, the signals result in the 
following lateral range function: 

p(x) = 1 - exp{-[2/(wr) ] • (r' - x^)^) for |x| < r 

(18) 

p(x) = 0 for |x| > r. 




Figure 3. The encounter geometry for the two intermittent signal 
models described here. 
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For signals whose occurrence is determined by a Poisson 
process and for which 5t > 0, signals can overlap. If this is 
allowed, then Equation 18 can be modified to describe this case 
by adding 6t/r to the term in the exponent of Equation 18 that 
is within the square brackets. In particular, note that this 
modified Equation 18 can be approximated by the bottom equality 
in Equation 17 when (2/w»r)(r* - x^ ) ^ + St/r « 1. This 
implies that when the expected time r between signals is large 
relative to the exposure time (2/w) (r^ - x^ ) ^ + 6t, the periodic 
occurrence model and the random occurrence model are effectively 
equivalent. 
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X. A RandoB Search Model 



A search of a region in which a target moves on a track 
that consists of a number of straight segments placed in such a 
way that in a limiting sense every section of the region is 
equally likely to be searched on a segment is referred to as a 
random search in Reference 6. Representation of a search region 
with the track segments that could be imagined to be the tracks 
of a random search are shown in Figure 4 . 




Figure 4 . A search region and a track that could be described as 
a random search track. 

Two developments of a model to describe this kind of search 
are contained in this section. The first development is based on 
the following conditions; 1. A target is at a fixed position 
within a defined search region. 2. A searcher's track is a 
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sequence of straight line segments that are within the search 

region. 3. The searcher's detection system is such that while on 

a track segment, a rectangle is searched that is contained within 

the search region, is of length equal to the length of the track 

segment and is oriented so that its long axis is parallel to the 

track segment. 4. The probability that the searcher's detection 

system will detect a target while on a track segment with a 

search rectangle that does not contain the target is zero. The 

probability that the searcher's detection system will detect a 

target while on a track segment with a search rectangle that 

contains the target is p(x) where x is the target's lateral 

range for the track segment and p(x) is the lateral range curve 

for a complete straight line encounter lateral range x. A 

representation of a search rectangle is shown in Figure 5 below. 

5. The track segments are located in such a way that the event 

that the target is within the search rectangle associated with a 

track segment is independent of the event that the target in the 

search rectangle associated with any other track segment. And 

the probability of the event is equal to the ratio of the area of 

<• 

the search rectangle to the area of the search region and, given 
a target is within a search rectangle, its position is uniformly 
distributed over the rectangle. 

Condition 4 implies that the random search model is based 
on the concept of a complete straight line encounter. The 
definition of an encounter that is intended here is that given in 
Section VI. This implies that in the random search model the 
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time to resolve a false alarm is zero. However, for the model, 

P(j and pf are considered to be determined by some criterion 
such that Pf is less than one. Consequently, although the time 
to resolve a false alarm is zero in the model, the cost 
associated with a false alarm is not zero. (A concise model that 
accounts for the time to resolve false alarms is described in 
Reference 13.) Condition 4 also implies that when a searcher is 
on a track segment with a search rectangle that contains a 
target, the encounter is a complete straight line encounter. And 
Condition 5, which can be considered to specify a random 
arrangement of the track segments, implies that when this is the 
case, for the complete straight line encounter, the target's 
lateral range is a random variable that is uniformly distributed 
between -b/2 and b/2 where b is width of the search 
rectangle (the dimension of the rectangle perpendicular to the 
associated track segment) . 



Swept Area 



Track 

Figure 5. A track segment and its associated search rectangle, 
that could correspond to a search with an aircraft mounted 
infrared detection system. 
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Based on the above considerations, the probability that a 
target will be detected while a searcher is on a track segment 
with an associated search rectangle that contains the target is 
given by: 



where for -b/2 < x < b/2 and fx(^) “ 

p(x) = 0 otherwise. Note that the left side of Equation 19 
applies to any complete straight line encounter in which the 
target's lateral range for the encounter is considered to be a 
random variable with a distribution determined by the probability 
density function fx(^) • If it is not given that the target is 
within the search rectangle associated with a track segment, then 
the unconditional probability that the target will be detected on 
the track segment is given by: (W/b)*(6A/A) where 5A is the 

area of the search rectangle associated with the track segment 
and A is the area of the search region. With 1 the length of 
the rectangle, 6k = b»l and the probability becomes: (W*l)/A. 

Then, since the event that the target will be in the search 
rectangle of a track segment is independent of the event that it 
will be in the search rectangle of any other track segment, the 
probability p that a random search consisting of m track 
segments will detect the target is given by: 

1 - [1 - (W-li)/A][l - (W«l 2 )/A] ••• [1 - (W*ln)/A] where li 
is the length of the i^^ track segment. The probability is 
also given by: p = i - exp{S ln[l - (W«li)/A]} where the sum 



(19) 
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index i = 1,2, ••• ,n. If (W*li)/A « 1 for i = 1,2, ••• ,n, 
then this expression can be approximated by: 



where 1 = E is the track length of the search. Equation 20 

is known as the random search formula. 

The second development of the random search formula is based 
on Equation 5 and a detection rate for a random search given by: 
T(t) = W*v(t)/A. With this detection rate and Equation 5, the 
random search formula is given by: 

(20a) P(T < t) = 1 - exp {-[W*l(t) ]/A) 

where l(t) is the track length for a random search that starts 
at time 0 and ends at time t and 



Replacing P(T < t) by p and l(t) by 1 gives Equation 20. 
In the form of Equation 20a, the random search formula indicates 
explicitly the relation between the probability of detection and 
the duration of a random search. Note that Equation 20a implies 
that the sweep width is independent of speed over the range of 
speeds in the encounter. 

Reference 14 contains an example of an application of the 
method used in the second development of the random search 
formula to a random search where the search region expands with 
time. 



( 20 ) 



p = 1 - exp[-(W*l)/A] 



(20b) 



l(t) 
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XI. Ladder and Barrier Search Models 

In some barrier searches, the barrier search track is a 
ladder search track relative to a reference system that moves 
with the target. This fact is used in the barrier search model 
development that follows the two ladder search model developments 
below. The first ladder search model is referred to as an ideal 
ladder search model because of the idealizations that are 
involved in its description of a ladder search. The second 
ladder search model is referred to as a degraded ladder search. 

It can be considered to describe a ladder search track in which 
navigational errors result in omissions and overlaps in coverage. 

An Ideal Ladder Search Model: The model is based on the 

following conditions; 1. A ladder search region is a rectangle 
that contains a fixed target. 2. During a search of the region, 
the searcher's detection system moves on a set of m parallel 
track segments of length b separated by a distance s. 3. As 
the detection system moves along a track segment, it searches a 
rectangular strip of length b and width s within the search 
region. 4. The m rectangular strips that correspond to the m 
track segments completely cover the ladder search region with no 
overlap. 5. If a target is within the rectangular strip 
corresponding to a track segment, then there will be a complete 
straight line encounter between the target and the detection 
system when the detection system moves along the track segment 
and the lateral range of the encounter will be uniformly 
distributed across the width of the strip. If the target is not 
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in the rectangular strip, then there will not be an encounter and 
the probability that the target will be detected while the 
detection system is on the track segment is zero. 




Figure 6. A schematic representation of a ladder search geometry 
for a case in which the ladder search track segments are 
superimposed on and bisect their corresponding rectangular 
strips. 

Since targets outside of the rectangular strip that 
corresponds to a track segment cannot be detected while a 
detection system is on the track segment because of Condition 5, 
in the model, the sweep width W of a searcher's detection 
system must satisfy the relation W < s . In particular, W = s 
only holds when the detection system detects a target that is in 
a rectangular strip with probability one for any target lateral 
range. This kind of detection system is sometimes referred to as 
a cookie cutter detection system. However, this terminology can 
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be misleading since it suggests the detection system detects 
equally well for all azimuths. But this is not a requirement on 
the system in order that W = s. 

The ideal ladder search model implies that if the conditions 
of the model are satisfied, then the probability p that a 
target will be detected by a an ideal ladder search is given by: 

(21) p = W/s 

where W/s < 1. The quantity W/s is called the coverage 
factor in this case. 

A Degraded Ladder Search Model: The above model implies 

perfect navigation in addition to other idealizations. A model 
of a ladder search is given in Reference 6 that could be used 
for cases in which this is a poor assumption. The model which is 
referred to here as a degraded ladder search model can be 
considered to describe navigational inaccuracies in terms of 
omissions and overlaps of the rectangular strips. It can be 
developed as follows: Consider a random search in the ladder 

search region whose track length is equal to the search track 
length required to complete an ideal ladder search, that is, a 
track length 1 = m*b. The degraded ladder search model 
describes the result of omissions and overlaps in a ladder search 
to be such that the probability of detection for this random 
search is equal to the probability of detection for the degraded 
ladder search. Consequently, since the area of the ladder search 
region is m*s*b, for the degraded ladder search model: 

(22) p = 1 - exp(-W/s). 
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Here, the requirement that the coverage factor W/s < 1 for 
Equation 21 can be relaxed. However, it should still be 
considered as an approximate condition. 

The condition that the target be fixed within the 
rectangular search region is critical to both Equation 21 and 
Equation 22. However, these results are also applicable to a 
search for a moving target under the conditions that are 
described next. 

A Barrier Search Model: A target moves with a constant 

course and a constant speed u. Both the target's course and the 
target's speed are known by a searcher. The searcher establishes 
a barrier of width b that is perpendicular to the target's 
track and moves on the barrier with a speed v > u. The barrier 
is designed so that in a reference system relative to the target 
the barrier search is a ladder search that satisfies the 
conditions for a ladder search that are given above. There are 
two cases to consider: 1. The barrier is established in front of 
the target. 2. The barrier is established behind the target. 

From the search geometry for a barrier established in front 
of the target, it can be seen from Figure 7 below that 
6 = sin”^(u/v) and d = vt where t = s/(v + u) is the time to 

move from one search leg to the next. The angle 0 and the 

perpendicular distance d that depend on u, v and s, and the 

width of the barrier b are the quantities that are required in 
order to establish the barrier operationally. 



43 



Target 
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Figure 7. A barrier search track shown for a barrier established 
in front of the target. The track is shown in a reference system 
fixed relative to the earth. 

For a barrier that is established in front of a target, one 
of three barrier types will result. A barrier's type is 
determined by the relation of the distance d to the distance 

g = ut where the time t - h/ (v' - u' ) ^ is the time to complete 

a search leg (cross the barrier) . The barrier type is determined 
as follows: 1. For g < d, the barrier is an advancing barrier. 

2. For g = d, the barrier is a stationary barrier. 

3. For g > d, the barrier is a retreating barrier. 

For a barrier established behind the target, there is only 
one barrier type and it is called an overtaking barrier. For an 
overtaking barrier, 0 = sin“^(u/v) as for a barrier established 
in front of the target. But, for an overtaking barrier, 

T = s/(v - u) and d = vs/(v - u) . 
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Given the target crosses the barrier, the probability of 
detection for an ideal barrier search is given by Equation 21 
and the probability for a degraded barrier search is given by 
Equation 22 where the terminology refers to the nature of the 
ladder search in the reference system moving with the target. 

A discussion of an application of these two equations to a search 
for a magnetic anomaly target is given in Reference 15. 
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XII. A Target State Estimation Procedure 

A target state estimation procedure based on bearing 
observations is developed in this section that generates point 
estimates of a target's position and velocity vector coordinates 
in a rectangular coordinate system. The procedure is based on a 
model in which bearing errors are unknown and are not determined 
by random variables with known distributions. Because of this, 
confidence regions for the estimates are not generated by the 
procedure. However, for a moving target, it illustrates general 
characteristics of bearings only target motion analysis (THA) . 

The model is defined as follows: 1. The target moves in a plane 

with a constant but unknown course and speed. 2. Observations of 
the target are made from known positions at known times. 3. The 
observations provide only target bearings with unknown errors. 

The model geometry is shown in Figure 8 . 
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Figure 8. The geometry of the target motion analysis model. 
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The procedure criterion is: For observations from n positions, 

chose target position estimate and target velocity component 
estimates Uj^ and Uy that make the sum S of the squares of 
the algebraic distance between the estimated position and the 
bearing line corresponding to the observed bearing a minimum. 

From Figure 8, the algebraic distance can be written as 
di = [x(i) - Xo(i)]*cos - [yt(i) - yo(i)]*sin Because of 

the requirement that the target move with constant course and 
speed during the encounter, the number of independent estimates 
is reduced from 2n to 4: Xt(l) , Yt(l)» 'Jx These 

four estimates define a target's motion by the relations: 
xt(i) = xt(l) + Ux* (ti - ti) and yt(i) = Yt(l) + ^y* (^i ~ t^) . 
To determine estimates of the target state parameters, take the 
partial derivative of the sum S = E (dj^)’ with respect to each 
of them. Then set the four partial derivatives equal to zero. 
This creates four linear equations in x^Cl) / yt(^)» *Jx ^y 

whose solution are the desired estimates xt(l), '^x 

Uy. In matrix notation, the equations can be represented by: 

AX = B where the elements of X are x^(l), yt(i)» ’^x' 
uy^ A necessary condition for a unique solution for X is that 
n > 4. Otherwise, the determinant of A will be equal to zero. 
The procedure can also be used if a target's course and speed are 
constant and known and, in particular, if the target is 
stationary so that u^ and Uy both equal zero. In this case, 
since the number of unknowns is two, the number of linear 
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equations is also two and a necessary condition for a unique 
solution is n > 2. 

Now, suppose the observations are at positions and times 
that correspond to the positions and times of an observer moving 
on some constant course at some constant speed (including zero 
speed) . In this case, the observation position coordinates are 
related by the following equations: XQ(i) = Xo(l) + v^Ct^ - t^) 

and yo(i) = yo(l) Vy(ti - ti) where v^ and Vy are the 
required velocity components of the observer. Using these 
equations of motion, the matrix equation AX = B can be 
transformed to the matrix equation AX'= 0 where the elements 
of the matrix X are related to the elements of the matrix X' 
by the equations: " ^qCI)/ ^^21 = YtC^) ” YoC^)/ 

*31 = ^x ~ Vx and X41 = dy - Vy. 

Since the linear equations represented by AX' = 0 are 
homogenous, they do not have unique solutions and consequently 
neither do the equations represented by AX = B. However, if 
there is at least one observation whose time and position is not 
determined by the above equations of motion, then the 
transformation from X to X' cannot be made, and in general a 
unique solution for X can be found. If the observations are 
made from a platform that is moving with a constant course and 
speed, this condition can be achieved by either changing the 
course, the speed or both prior to completing the observations. 

Estimation models that describe bearing error as a random 
variable provide a basis for determining confidence regions for 
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point estimates. A model is developed in Reference 16 that does 
this for either target bearing observations made from two or more 
points simultaneously or for a target that is stationary relative 
to the observation points. 
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XIII. Position Distributions That Change with Motion 

In this section, target motion models are discussed 
that provide a basis for determining position distributions that 
change with time. In the first model, a target moves in a plane 
with a constant course and speed. In the second model, a target 
moves in a plane with a constant course and speed during a time 
step and at the end of the time step changes either its course or 
speed or both. 

Position Distributions based on the First Motion Model: Two 

cases of this model are considered that provide the basis for 
determining a position distribution analytically. In both cases, 
the target's position coordinates in a rectangular coordinate 
system are determined by a circular normal distribution at time 
zero. And, in both cases, during the motion the target's 
velocity does not change. 

The First Case: In this case, the targets 's coordinates 

X(t) and Y(t) at time t are determined by the vector sum 
with components: X(t) = X(0) + Ujj»t and Y(t) = Y(0) + Uy*t 

where X(0) and Y(0) are the target's coordinates at time zero 
and Ujf and Uy are the target's velocity components. The 
coordinates X(0) and Y(0) are determined by a circular normal 
distribution whose mean vector is (0,0) and whose standard 
deviation is a. And the velocity components Ujf and Uy are 
determined by a circular normal distribution whose mean vector is 
(ttjji'iy) ^nd whose standard deviation is cTy. The random vectors 
[X(0),Y(0)] and (Ux*t,Uy»t) are independent and both are 
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normal random vectors. Since they are normal random vectors, the 
random vector [X(t),Y(t)] is also a normal random vector. And, 
since they are independent, its characteristic function is equal 
to the product of their characteristic functions. This implies 
that [X(t),Y(t)] has a circular normal distribution whose mean 
vector is (ttx*t,Uy*t) and whose variance is (a^ + o{j*t^ ) . For 
a more complete discussion of the basis for the above argument, 
see Reference 7. 

The Second Case: The procedure that is used to determine 

the distribution of the random vector [X(t),Y(t)] for the 
second case can also be used to determine it for the first case. 
It is based on the following consideration: In general, the 

joint density function of a target's coordinates at some time t 
greater than a zero of time is determined as follows: 



(23) fx(t) ,Y(t) = -oofx(O) ,Y(0) (q/s;0) fv,w(v,w) dvdw 

where q = x - vt and s = y - wt and V = Uj^, W = Uy, v = u^ 
and w = Uy. Equation 23 can be developed as follows: To first 
order, f^^t) ,Y(t) probability that a 

target's coordinates are in an element of area <5x*<5y and for 
given values of v, w and t, target positions in an element of 
area <Sq*5s will be translated to an element of area 5x*5y 
that is identical in form and size to 6q*<5s. And so, to first 
order, fx (O) » Y (0) ® ' °) * *^^* ^y * ^V, is 

probability that the target's coordinates at time 0 are in an 
element of area (5q*6s that is located such that the target's 
coordinates will be in the element of area 5x*5y at time t 
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since x = q + vt and y = s + w*t. And, to first order, the 

sum of such probabilities for all pairs of values of v and w 
is also the probability that the target's coordinates at time t 
are in the element of area Sx’Sy. In the limit after equating 
the two expressions for this probability and cancelling the 
common factor 6x*Sy, Equation 23 results. 

In the first case, the joint distribution of the velocity 
components (U^/Uy) is circular normal with mean vector 
(Ujj,Uy) and standard deviation and the joint distribution 

of [X(0),Y(0)] is circular normal with mean vector (0,0) and 
standard deviation a. In the second case, the target's speed 

u is known but its course is equally likely to have any value 

between 0 and 2ir , so Ujf = U‘sin § and Uy = U‘Cos * where ♦ 

is a random variable with a uniform distribution and values 
0 < 0 < 27T. In this case, only the single random variable * is 
required to determine the target's velocity. With the velocity 
components in polar coordinates, the integral of Equation 23 is 
a single integral over <p and the integrand of the integral is 
(l/27ra* ) •exp[-(q* + s^ )/2a' ] • (l/2w) where now q = x - u*t*sin <p 
and s = y - u*t*cos <p. The result of the integration is: 

(24) (l/27ra^) exp{-[x^ + y^ + (u»t)^]/2aM + yM^*u*t/a^] 

where Iq indicates the hyperbolic Bessel function of zeroth 
order. Expression 24 is plotted in Reference 6 for several 
values of t. There, a single variable r is substituted for 
(x* + y^)^, the the target's range from the origin. The plot of 
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the expression shows a limiting characteristic of the 
distribution that can be indicated as follows: First replace 

(x* + y* ) "* by r in Expression 24. Then, by first multiplying 
and then dividing the resulting expression by exp(-r*u*t/a ^ ) • 
the following expression can be obtained: 

(25) 1/ (27ra^ ) exp{ -[1/ (2a^ ) ] (r - u* t) ^ ) Iq (r ‘U't/a^ ) exp(-r *u*t/a^ ) . 
The second factor in Expression 25 is proportional to the 
density function of a normal random variable whose mean is u*t 
and whose standard deviation is a. And the last two factors are 
such that their product changes slowly as a function of r 
relative to the second factor. The consequence of this is that a 
plot of Expression 25 against r for values of t greater than 
4*a/u has the appearance of a normal density function. In 
particular, this indicates that u*t in the exponent of the 
second factor in Expression 25 can be viewed as the radius of an 
average furthest on circle. 

A target's rectangular coordinates and its range and bearing 
from the origin are related by: X = R sin 9 and Y = R cos 9. 

By using these equations. Expression 24 can be transformed to: 

(26) (1/2W) (r/27raMexp{“[r* + (ut) * ]/2aM I q (^ut/a* ) 

which is the joint density function of the random variables R 
and 9. The marginal density function of R can be obtained 
from Expression 26 by integrating this joint density function 
over the possible values of 9 which in radians is over the 
interval 0 to 2n . Consequently, the marginal density for R 
can be obtained from Expression 26 by multiplying it by 2n. 
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Tabulated values of the cumulative distribution function Fp(r;t) 
for the marginal distribution of R for the second case are 
listed in Reference 17. 

Position Distributions Based on the Second Motion Model: In 

the second motion model a target's motion is not restricted to a 
constant course and speed throughout the motion. At the end of a 
time step, a target's course or speed or both can change. In 
general, this leads to a position distribution that can not be 
described analytically. And, in general, a monte carlo 
simulation method is required in order to determine a position 
distribution that is based on such a model. As an example of a 
case in which the distribution can be described analytically, see 
Reference 18. As an example of a case in which it can not, 
suppose a target's initial position is described in terms of a 
number assigned to a subregion in the xy-plane where the number 
assigned represents the probability that the subregion contains 
the target at an initial time. In addition, suppose for each 
subregion a course and speed distribution is determined by 
assigning numbers to course and speed pairs where a number 
represents the probability the target will have the course and 
speed at the initial time given it is in the subregion at that 
time. Next suppose for each course and speed pair there is a 
time distribution that determines the duration pf the course and 
speed pair and that the time distribution is determined by a 
number assigned to each discrete time point where the number 
represents the probability that the target's course and speed 
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pair will be determined by a new course and speed distribution. 

By extending this kind of procedure and then implementing it in a 
monte carlo simulation, one can generate complex position 
distributions that describe a target's position at discrete time 
points in terms of probabilities assigned to subregions of a 
region in which the target is contained. 
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XIV. Position Distributions That Change with Search 

Suppose a target’s position at some time is described 
by a position distribution. Now suppose information becomes 
available that a search has been conducted for the target and 
that the target has not been detected. Or suppose the 
information is that the target has been detected. In the first 
case, negative information is available that can be used to 
modify the position distribution. In the second case, positive 
information is available that can be used to modify the position 
distribution. In both of the cases, a target's position 
distribution is assumed to be specified in terms of a set of 
discrete probabilities where each probability corresponds to a 
subregion of the region that contains the target and each is the 
probability that the target is in that subregion. 

Position Distributions and Negative Information: For a 

region that contains a target and consists of n subregions, let 
the event Sj^ = {the target is in the i^^ subregion). And let 
the event C = (no contact). Then, given no contact in a search 
of the region, the targets 's position distribution can be 
modified as follows: 

(27) P(Si|C) = P(c|Si) -P(Si)/P(C) 

where i = 1,2, ••• ,n and P(C) = Z P(C| Sj ) *P(Sj ) with the sum 
from j=l,2, ••• ,n. Note that Equation 27 can be obtained by 
using Bayes theorem. To illustrate how Equation 26 might be 
used, suppose that a search in a subregion is considered to be a 
random search and that the sweep width of detection system 
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against the target depends on subregion being searched. For this 
case, let be the area of the i^^ subregion and let be 

the sweep width in that subregion. Then, given no contact in a 
search of a subregion, P(C|Si) = exp[-(Wi* li)/Ai] where 1^ is 
the track length of the searcher in the i^^ subregion. Given 
values for P(Si) , <5^ and 1^ for i=l,2, ••• ,n a position 

distribution can be determined that has been modified by the 
negative information. 

Position Distributions and Positive Information: In the 

case of positive information, the event C = {a contact) occurs. 
This event is the union of two mutually exclusive and exhaustive 
events: T^ = {a true contact) and = {a false contact). 

Relative to the Venn diagram of Figure 1, the event Sj^ defined 
above corresponds to C corresponds to T^, 

corresponds to (D^ n and F^, corresponds to (Dj^ n Hq) . 

For generality, suppose true contacts do not localize a target to 
a single subregion. Then after a search of a region that has 
resulted in a contact, the target's position distribution can be 
modified as follows: 

(28) P(Si|C) = P(Si|Tc) *P(Tc|C) + P (Si | Fc) • P ( Fc I C) 

where i = 1,2, ••• ,n since T^ = T^ n C, F^ = F^ d C and 

P(Si|C) = (P[(Si n Tc) + P(Si n Fc)])/P(C). The probability 

p = P(Tc|C) has been called the credibility of the contact. In 
terms of p Equation 28 becomes: 

(29) P(Si|C) = P(Si|Tc)-p + P(Si|Fc)-(l - P) . 
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For the model, P(Si|Tc) can be determined by the coverage 
characteristics of the detection system that is used to make the 
contact. In particular, P(Si|Tc) = P(Tcl S^) *P(Si)/P(Tc) where 
P(Tc) = 2 P(T^ I Sj ) *P(Sj ) and the sum index j = 1,2, ••• ,n. 
Then, using the correspondence between the events Sj^ and 
and the events and n , a definition for P(Tc|Sj^) 

is P(Tc|Si) = (Pd)i where (Pd)i is an average probability of 
detection over the i^^ subregion. In one positive information 
model, P(Si|F) = P(Si) (false contacts supply no information 
about a target's location), and p is determined subjectively 
based on factors associated with the detection system used to 
make the contact. The probability determined by 1 - p has been 
called the false alarm probability. However, it is not the 
probability pf = P(Dj^|Ho). Since 1 - p = P(F<^|c), it is the 
probability P(Ho|Di). 

The above approach can also be used to develope a procedure 
to modify a position distribution given negative information. To 
do this, note that the event C is the union of two mutually 
exclusive and exhaustive events: T^ = {a true no contact) and 

F^ = {a false no contact). Relative to the Venn diagram of 
Figure 1, C corresponds to Dq, T^ corresponds to (Dq n Hq) 
and F^ corresponds to (Dq n H 2 ^) . In analogy to Ecjuation 29, 
(30) P(Si|C) = P(Si|Tc) -P + P(Si|Fc) • (1 - P) 

where p = P(Tg|c). Relative to the Venn diagram of Figure 1, 

P = P(Dq|Hq). In analogy to the choice for P(Sj^|T<;.) above, one 
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could choose P(Si|Tc) = P(Si) and P(Si|Fc) could be 
determined by the coverage characteristics of the detection 
system. In particular, P(Si|Fc) = PCF^I Sj^) •P(Si)/Fc) where 
P(Fc) “ ^ P(Pcl ) *P(Sj ) index j = 1,2, ,n. 

Using the correspondence between the events Fq and Dq n 
one can write P(Fc|Si) = 1 - (Pd)i. 

Both in the case of this procedure for negative information 
and the procedure for positive information, (Pcj) i can be viewed 
as a measure of the specified performance of a detection system 
and both p and p can be viewed as measures of the degradation 
in its specified performance do to various operational factors. 

In particular, for p = 0, the expression for P(Sj^|C) is the 
same for the two negative information procedures if P(C|Sj[) is 
determined in the same way for both cases. In the view just 
expressed, this implies that the first procedure applies to a 
search in which a detection system is operated so that it 
achieves its specified performance. In the case of the positive 
information procedure, if p = 1, then the expression for P(Si|C) 
is the same as that which can be obtained from 
(31) P(Si|C) = P(C|Si) *P(Si)/P(C) 

where P(C) = S P(C | Sj ) ♦ P(Sj ) and the sum index j = 1,2, ••• ,n 
if P(C|Sj^) is determined in the same way for both cases. 
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XV. Search Models and Search Theory 

Search theory provides a basis for determining optimal 
search plans for a target whose motion and location are 
determined within some bounds. Here, an optimal search plan is 
one for which the probability of finding a target within a given 
length of time is a maximum, the expected time to find a target 
is a minimum given the target is found or a search plan for which 
some other optimal search criterion is satisfied. 

Search theory results are based on models of the search 
process. To the degree that a search model describes a search 
process, an optimal search plan for a target that is based on the 
search model should provide guidance for the development of an 
operationally feasible search plan. However, because of the 
limitations of analytical search models, an optimal search plan 
that is based on an analytical search model may give only initial 
guidance in this regard. The optimal search plans that are 
described below illustrate this. The search plans are based on 
the random search model. Because of this, the requirement on the 
location of search track segments is not realizable and the time 
to resolve false alarms is ignored. 

Optimal search plans based on search models implemented 
through a monte carlo simulation are not considered here. 

However, with sufficient information, such plans have the 
potential of being both implementable and more optimal in a real 
sense than an optimal search plan based on an analytical search 
model . 
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Three Optimal Search Plans: The three optimal search plans 

differ through their definition of optimality. However, each one 
is based on the following search model: A target is fixed at 

some point in a region that consists of n subregions. A search 
in a subregion is a random search in the sense of the definition 
in Section X and a searchers sweep width there is a constant. 

In addition, a search of a subregion will not detect a target 
which is in another subregion. To determine a plan, let 
S^ = {the target is in subregion i} for i = 1,2, ••• ,n and 
let pj^ = P(Si) be the prior probability that the target is in 
the i^^ subregion. Let be the sweep width in the i^^ 

subregion. Let ^i = A^/Wj^ where Aj^ is the area of the i^^ 

subregion and is the expected track length to find the 

target by a search of the i^^ subregion given the target is in 
the i^^ subregion, a characteristic length. The probability P 
that the target will be detected by a random search is given by: 

(32) P = 2 [1 - exp(-li/<Si) ] -Pi 

where the sum index i = 1,2, ••• ,n and Ij^ is the track 

length of the search in the i^^ subregion. 

The first optimality criterion is: Choose Ij^ so that P 

is a maximum subject to the two constraints: 1. 1 = S 1^ and 
2. l£ > 0 where the index i = 1,2, ••• ,n. Determining this 
choice is a nonlinear optimization problem whose solution is 
given in Reference 19. It is: 

l+/<Si = ln(pi/6i) - L(k) i = 1,2, ••• , k 

(33) 

li/5i =0 i = k+l,k+2, ,n 
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where L(k) = ( 1/E<S j ) •!:[ 5j • In (pj/<S j ) ] + l/r«Sj and the sum index 
j = 1,2, ••• ,k, where the subregions are relabeled so that the 
following order relation holds: Pi/<Si > P2/^2 > *** > Pn/*^n 

where k is chosen so that for k+1 the solution for 1)^+1 
using L(k+1) is either negative or zero. 

The second optimality criterion is: Choose Ij^ so that P 

is a maximum subject to the two constraints: 1. c = S Cj^ and 

2. c^ > 0 where the index i = 1,2, ••• ,n, cj^ = k^'l^ is the 
cost of the search in the i^^ subregion and kj^ is the cost per 
unit track length in that subregion. For this criterion, the 
solution to the corresponding nonlinear optimization problem can 
be obtained from Equation 33 by replacing by = ^^i*^i 

and labeling the subregions so that Pi/^i > P2/^i > ••• > Pn/^i* 
The basis for this can be seen by replacing li/^i by its 
equivalent in the exponential term in Equation 32. 

The third optimality criterion is: Choose 1^^ so that the 

expected utility of finding the target is a maximum subject to 
the two constraints: 1. 1 = S Ij^ and 2. Ij^ > 0 where the 
index i=l,2, ••• ,n. For this criterion, the solution to the 
corresponding nonlinear optimization problem can be obtained from 
Equation 31 by replacing pj^ by q^ where q^ = u^'p^ and uj^ 
is the utility of finding the target given it is in the i^^ 
subregion. And, in addition, labeling the subregions so that 
qi/5i > q2/^2 > *•* > qn/*^n* basis for this can be seen by 

multiplying the summation term in Equation 32 by uj^ so that 
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the resulting equation gives the expected utility of finding the 
target with the search. 

Equation 33 can be used to determine an order of search for 
the subregions which will effectively minimize the expected track 
length required to detect a target given it is detected. To do 
this, divide the available track length 1 into units small 
enough so that with a single unit only the 1®^ subregion would be 
searched. Then allocate one unit to the search of the 1®^ 
subregion. If the search is unsuccessful, determine the optimum 
allocation for two units. Then search with a second unit so that 
the first search with the first unit plus the second search with 
the second unit satisfy the optimum allocation for two units. If 
the search is unsuccessful, continue in this fashion until either 
the target is found or all the track length is expended. That 
this allocation order will effectively minimize the expected 
track length required to detect a target given it is detected can 
be argued as follows: Let L be the track length at detection, 

let 1^ be a unit of track length and let n be the number of 
units. Then the value of the probability P(L < i*lu) that the 
target will be detected on or before the i^^ step of the search 
for the given allocation order will be greater than or equal to 
its value for any other allocation order with the same allocation 
step size. Since the value of P(L < 1) will be equal to its 
value for any other allocation order of the optimum allocation 
and since P(L < i*ly|L < 1) = P(L < i*l,j)/P(L < 1), the value of 
the distribution function FL(i*lu|L < 1) = P(L < i'l^lL < 1) 
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will be greater than or equal to its value for any other 
allocation order. This implies that the expected track length 
given detection E(l|l < 1) = E [1 - FL(i*lu|L < 1)] where the 
sum index i=l,2, ••• ,n is effectively a minimum for the 
given allocation order. A search based on the optimum allocation 
given by Equation 33 and the given allocation order is 
equivalent to the following search: After an allocation of track 

length l^ and an unsuccessful search, new values for P(Si) are 
calculated using Equation 27 and then Equation 33 is used with 
these new values to determine the next optimum allocation. A 
discussion of this procedure is given in Reference 6. And an 
example of its application is given in Reference 20. 

Equation 33 also defines an optimal search plan for a 
detection system that searches beams and can be described by 
Equation 31 by replacing 1^ by tj^ where t^ is the time the 
i^h beam is searched and by replacing (5^ by where r^, a 

characteristic time, is the expected time to detect the target by 
a search of the i^^ beam given the target is in the i^^ beam. 

For a more extensive discussion of search theory and its 
application to military operations research, see Reference 21. 
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